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ABSTRACT 



Aims. A multi-transition survey of HCN (sub-) millimeter line emission from a large sample of AGB stars of different chemical 
type is presented. The data are analysed and circumstellar HCN abundances are estimated. The sample stars span a large range of 
properties such as mass-loss rate and photospheric C/O-ratio. The analysis of the new data allows for more accurate estimates of the 
circumstellar HCN abundances and puts new constraints on chemical models. 

Methods. In order to constrain the circumstellar HCN abundance distribution a detailed non-LTE excitation analysis, based on the 
Monte Carlo method, is performed. Effects of line overlaps and radiative excitation from dust grains are included. 
Results. The median values for the derived abundances of HCN (with respect to H 2 ) are 3 x 10~ 5 , 7 x 10~ 7 and 10~ 7 for carbon stars 
(25 stars), S-type AGB stars (19 stars) and M-type AGB stars (25 stars), respectively. The estimated sizes of the HCN envelopes are 
similar to those obtained in the case of SiO for the same sample of sources and agree well with previous results from interferometric 
observations, when these are available. 

Conclusions. We find that there is a clear dependence of the derived circumstellar HCN abundance on the C/O-ratio of the star, 
in that carbon stars have about two orders of magnitude higher abundances than M-type AGB stars, on average. The derived HCN 
abundances of the S-type AGB stars have a larger spread and typically fall in between those of the two other types, however, slightly 
closer to the values for the M-type AGB stars. For the M-type stars, the estimated abundances are much higher than what would 
be expected if HCN is formed in thermal equilibrium. However, the results are also in contrast to predictions from recent non-LTE 
chemical models, where very little difference is expected in the HCN abundances between the various types of AGB stars. 



Key words. Stars: AGB and post-AGB 
abundances 



1. Introduction 



Stars: carbon - Stars: late-type - Stars: mass-loss - (Stars:) circumstellar matter - Stars: 



The final stages of the evolution of low- to intermediate-mass 
stars are, to a large extent, dominated by considerable mass loss 
in the form of an intense stellar wind. The circumstellar enve- 
lope (CSE), created by the wind as the star ascends the asymp- 
totic giant branch (AGB), contains the products of the elements 
produced by internal nuclear processes, and dust grains formed 
in the inner part of the wind near the photosphere. Hence AGB 
mass loss contributes to the chemical evolution of the cosmic 
gas. In the outer CSE most molecules are photo-dissociated by 
the interstellar UV radiation field. The dust grains survive the 



* This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX) telescope, the IRAM 30 m tele- 
scope, the James Clerk Maxwell Telescope (JCMT), the Swedish- 
ESO Submillimeter Telescope (SEST), and the Onsala 20 m tele- 
scope. APEX is a collaboration between the Max-Planck-Institut fur 
Radioastronomie, the European Southern Observatory (ESO), and 
the Swedish National Facility for Radio Astronomy, Onsala Space 
Observatory (OSO). IRAM is supported by INSU/CNRS (France), 
MPG (Germany) and IGN (Spain). The JCMT is operated by the Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. The 
Onsala 20 m telescope is operated by OSO. The SEST was operated 
jointly by ESO and OSO. 



radiation and there is the possibility that some of the more re- 
fractory molecules, such as e.g. SiO and SiS, are incorporated 
into grains and protected from photodissociation. 

It is important to establish an understanding of the circum- 
stellar molecular line emission from AGB stars and to derive ac- 
curate physical and chemical properties of the CSEs. The mass 
return of such stars, and hence their contribution to the galac- 
tic chemical evolution, is dominated by the rare high-mass-loss- 
rate objects. These are highly obscured at shorter wavelengths by 
their large amounts of circumstellar dust and their photospheric 
abundances can not be determined using classical methods, such 
as visual and near-IR spectroscopy. Instead, one must rely upon 
estimates based on circumstellar emission. 

The structure of CSEs is relatively poorly known both 
on larger and on smaller scales. Single-dish observations do 
not generally resolve the sources and high resolution inter- 
ferometric observations are limited at present. Much of the 
current knowledge stems from interferometric observations 
of th e nearby, high-mass-los s -rate carbon star IRC + 10216 
(e.g.. iBieeing & Tafallal fl993t iDaval & Biegind fl993l 11995b 
iLucas et alJl995HFong et alj2006tlSchoier et al.l2006al l2607b), 
and th e over-all picture of the morphology of A GB-CSEs is 
crude (iNeri et alJll998t ICastro-Carrizo et aDl2010l) . The size of 
the emitting region is a crucial parameter for determining abun- 
dances, but theoretical estimates are uncertain, making interfer- 
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ometric observations important. Resolved images also test some 
of the basic assumptions of the radiative transfer model, e.g., the 
validity of spherical symmetry and smoothness of the wind. 

Observed abundances and spatial distributions of the cir- 
cumstellar molecular line emission will provide important tests 
of current chemical models of the stellar atmosphere and 
the CSE. The abundances of some common molecules such 
as HCN in carbon stars, and SiO in M-type AGB stars, 
can be reasonably well accounted for in chemical models 
where the molecules are formed in thermal equilibrium in, 
or near, the photosphere, and with little or no further pro- 
cessing until they are eventually photo-dissociated by the 
ambient UV-field. The detection of carbon-bearing species 
such as HCN, C N, CS and H 2 CO in envel o pes of M-type 
AGB stars (e.g. iDeguchi & Goldsmith! 119851: ILindqvist et ail 
1988; lOlofssonet al.1 (l 99 lULindqvistet al.1 1 1992b and oxveen- 



Table 1. Telescope data relevant for the new observations of 
HCN. 



bearing molecules like H 2 Q, H2CO, C3O and SiO in envelopes 
of ca r bon stars (e.g.. iKeadv & Ridgwavl 11993b iMelnick et al 
200 It iFord et al l 120041; i Hasegawa et all l2006t ISchoier et al 
2006bt iTenenbaum et al.1 l2006t iNeufeld et all 1201 lb points to 



the importance of non-equilibrium chemical processes. Models 
in which peri odic shock waves p r opaga t e in the extended 
photosphere dWillacv & C herchneff 1998; Duari eJiD [1991 
Agundez & Cernicharo 2006; Cherchneff 2006, 2012) appear to 
explain many of these observed abundances, however, a more 
statistical study is still lacking. Other molecular abundances may 
also be significantly affected by the formation of dust grains in 
the inner part of the circumstellar envelope (e.g., Agundez & 
Cernicharo 2006, SiO and H 2 0; Schoier et al. 2006b, SiO). 

HCN h as been detec t ed in AGB CSEs of all chemical 
types (e.g., iBieging et alJ 12000). and its formation is sensi- 
tive to shocks passing through the gas, making it a goo d 
probe of the importan t chemical processes (ICherchne fil 12006 ) . 
ILindqvist et al.1 (l2000b observed a sample of four carbon stars 
in the HCN J = 1 -> and CN N = 1 -> lines, using the IRAM 
Plateau de Bure interferometer, in an attempt to constrain their 
abundance distributions. From detailed numerical modelling 
they found a qualitative agreement with current theoretical mod- 
els of circumstellar chemistry when it comes to the size esti- 
mates (though, the statistics are poor and sources with inter- 
mediate mass-loss rates are missing). The spatial extents of the 
HCN and CN envelopes derived from the radiative transfer anal- 
ysis of the observational data, are systematically larger than 
predictions from photodissociation modelling, by about a fac- 
tor of two. Observations of circumstellar HCN line e mission 
from t he two M-type AGB stars TX Cam and IK Tau bv lMarvell 
(2005) clearly showed that HCN is also of photospheric origin 
in M-type AGB sta rs, which resolved a long-standing debate. 
iMuller et all (l2008b observed HCN ./ = 8 -» 7 emission, using 
the SMA, towards the M- type AGB star W Hya. Their data cor- 
roborates the findings of Marvel (2005) as to the origin of the 
HCN emission, and also probes gas still accelerating in the inner 
wind of the star. The recent analysis of the HCN Herschel/HIFI 
data (in combination with groundbased data) of the S-star^- Cyg 
(ISchoier et al.l20llh also suggest a more intricate HCN distribu- 
tion. 

In this paper we present new multi-transition observations of 
HCN line emission towards a sample of AGB stars with different 
properties, and circumstellar fractional HCN abundances (with 
respect to H 2 , hereafter referred to as abundance) derived using 
a detailed radiative transfer model. Together with existing data, 
these observations will significantly extend the observed range 
of parameter space, and provide a solid statistical test of current 
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chemical models in and near the photosphere of AGB stars, as 
well as further out in the circumstellar envelopes. 

2. Observations and data reduction 

New multi-transition HCN line observations (see Table [TJ were 
performed during 1995-2008 using the Onsala 20 m telescope 
(OSO), the Swedish-ESO 15 m submillimetre telescope (SEST), 
the JCMT 15 m telescope, the APEX 12 m telescope, the IRAM 
30 m telescope, and the NRAO 12 m telescope at Kitt Peak, 
USA. 

The OSO, SEST, JCMT, and IRAM observations were made 
in a dual beamswitch mode, where the source is alternately 
placed in the signal and the reference beam, using a beam throw 
of about 11' (OSO and SEST) or 2' (JCMT and IRAM). This 
method produces very flat baselines. At APEX and NRAO, the 
observations were carried out using a position-switching mode, 
with the reference position located +3' and +10' in azimuth, re- 
spectively. Regular pointing checks were made on SiO masers 
(OSO and SEST) and strong CO and continuum sources (APEX, 
JCMT, IRAM, and NRAO). Typically the pointing was found to 
be consistent with the pointing model within as 3" for all tele- 
scopes, except at the NRAO where it was found to be consistent 
within as 5". All receivers were single sideband receivers. 

The raw spectra are stored in 7^-scale and converted to 
main-beam brightness temperature using T m b-T*jT] m b. T* is 
the antenna temperature corrected for atmospheric attenuation 
using the chopper-wheel method, and r]^ is the main-beam 
efficiency. The adopted beam efficiencies, together with the 
FWHM:a of the main beams (0 m b), for all telescopes and fre- 
quencies are given in Table [1] The uncertainty in the absolute 
intensity scale is estimated to be about +20%. In Table [1] the 
energy of the upper level involved in the transition (£u P ) is also 
given, ranging from 4 K for the J -I level up to 43 K for the 
J = 4 level, illustrating the potential of these multi-transition ob- 
servations to probe a relatively large radial range of the CSE (see 
Sect. EXTJ ), 

The data were reduced by removing a low (typically first) or- 
der polynomial baseline and then binned (typically to a velocity 
resolution of about 1 kms~') in order to improve the signal-to- 
noise ratio, using XSQ- The observed spectra are presented in 
Figs IA.H4A.4l and velocity-integrated intensities are reported in 



1 XS is a package developed by P. Bergman to reduce 
and analyse single-dish spectra. It is publicly available from 
ftp : //yggdrasil . oso . Chalmers . se 
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Tables lATTI & IA.2l The intensities are given in main-beam bright- 
ness temperature scale (T m b). In some of the carbon-star spectra, 
the vi = 1, J = 3 — > 2, and J = 4 — > 3 lines can be seen at slightly 
lower frequency than the corresponding ground-state lines (es- 
pecially apparent in the spectra of e.g., Y CVn and CW Leo). 
The same vibrationally excited line appears to be present in 
the J = 3 — > 2 spectrum of the M-type star R Leo, however, it 
does not appear exactly at the correct frequency, and when re- 
observed it was not detected. 

In addition to the new data presen ted here we have als o 
used HCN lin e inten si ties reported by loTofsson et ail (119931). 
Olof sson et ail (119981). [Buiarrab al et all (11994 . Biegin g et all 
(20001). and lWoods et all (120031) . as well as data from the JCMT 
public archival All the data used in the excitation analysis are 
included in Tables |A~T1&|AT21 



3. Radiative transfer modelling 

3.1. General assumptions 

The CSEs are assumed to be spherically symmetric, produced 
by a constant mass-loss rate (M), and expanding at a con- 
stant velocity (v e ). The molecular line excitation analysis is 
performed using a detailed non-LTE radiative transfer code, 
ba sed on the Monte Carlo method and described in detail 
in ISchoier & OlofssonI d200lh . The code has been extensively 
tested and benchmarked ag ainst a wide variety of mol ecular line 
radiative transfer codes in Ivan Zadelhoff et alJ d2002l) . The lo- 
cal line width is assumed to be described by a Gaussian and is 
made up of a micro-turbulent component with a Doppler width 
of 1.0km s _1 and a thermal component that is directly calculated 
from the derived kinetic temperature structure (see Sect. I3.3l l. 

The analysis includes radiative excitation through vibra- 
tionally excited states of CO and HCN. Both the central star (ap- 
proximated by a blackbody), and thermal dust grains distributed 
in the CSE, can provide radiation fields capable of populating 
vibrationally excited states, and thereby affect the excitation in 
the ground vibrational state. The addition of a dust component 
in the Monte Carlo scheme is straightforward, as described in 
ISchoieret al.ld2002l) . 

The best-fit model is found by minimizing the total x 2 de- 
fined as 



-z 



(■Anod ^obs) 



(i) 



where / is the velocity-integrated line intensity and cr the uncer- 
tainty in the measured value. The cr-value is usually dominated 
by the calibration uncertainty assumed to be +20%. In general, 
the data quality is good and an uncertainty of 20% is a reasonable 
estimate. However, it is possible that the error is larger in some 
cases, and for very noisy lines we have added a higher uncer- 
tainty, thereby assigning a lower weight to the noisy lines in the 
modelling. The summation is done over N independent obser- 
vations. When modelling the dust continuum radiative transfer 
(Sect. I3.2l i. Eq. [T]is also used to find the best-fit model, but with 
the flux density, F^, instead of /. 

3.2. Dust continuum emission modelling 

The dust temperature structure and dust density profile are ob- 
tained from detailed radiati ve transfer mode lling using the pub- 
licly available code Dusty (llvezic & El itzur 1997). Here we are 



mainly interested in the potential effect that dust grains have on 
the excitation of molecules. A full treatment of dust radiative 
transfer coupled with a dynamical model for the wind, in or- 
der to independently measure physical properties of the CSE, 
such as e.g. th e mass-loss rate, is beyond the scope of this 
paper [see e.g. ISchoier et alJ (1201 ll) for an example of such 
modelling]. In the modelling, where the spectral energy dis- 
tribution (SED) provides the observational constraint, the dust 
optical depth specified at 10 /im, rio, the dust temperature at 
the inner radius of the dust envelope, T&(r\), and the effective 
stellar blackbody temperature, T*, are the adjustable parame- 
ters. The SED is typic ally constrained by JHKLM photomet- 
ric data [from 2MASS, iKerschbauml (1 19991) . and Kerschbaum 
priv. corn."!, IRAS fluxes, and in some cases sub-millimetre 
data (iGroenewegen et al]|1993l) . As mentioned above, the best- 
fit model is found by minimizing Eq. Q] 

The total luminosity, L*, for regularly pulsating 
stars, is obtained from t he per iod-luminosity relations of 
Groenewegen & Whitelockl d 19961 for carbon stars) and 



Whitelock et al.l (11994 for M- and S-type AGB stars). In the 



case of non-regularly pulsating stars, a generic luminosity of 
4000 Lq is adopted. The distance, D, is then obtained from 
the SED-fitting once the luminosity of the star is known. For 
the sources where reliable Hipparcos parallaxes exist, the 
corresponding distance has been adopted and the luminosity is 
then, in turn, obtained from the SED modelling. Amorp hous 
carbon dust grains with the optical constants given in ISuhl 
(2000j) are adopted for the carb on stars and astronomical 
silicates dJusttanont & Tielenslfl992l) for the M-type AGB stars. 
For the S-type AGB stars, either carbon or silicate grains are 
assu med based on the IR AS low-resolution spectra classifica- 
tion (IVolk & Cohenl [19891) . For simplicity, the dust grains are 
assumed to have the same radius, a^-Q. 1 /mi, and the same 
mass density, pd = 2.0 g cirT 3 (for carbon grains) and 3.0 g cm -3 
(for silicate grains). The corresponding dust opacities, k v , are 
calculated from the optical constants and the individual grain 
prope rties using standard Mie theory (Bo hren & Huf fman 
119831) . 

The results from the dust continuum emission modelling are 
listed in Tables[3]and|4] Th e resul ts for the S-type AGB stars are 
taken from lRamstedt et all (120091) . 

3.3. CO line modelling 

The physical properties of the circumstellar gas, such as the 
density, temperature, and kinematic structures, are derived 
from radiative transfer modelling of multi-transition millime- 
tre and sub-millimetre CO line observations. The CO data 



used in the analysis hav e been presented in Schoier & OlofssonI 
(120011 Olofss on et ail d2002l). bonzalez Delgado et alJ d2003|) . 



2009). 



http : //www . j ach. hawaii . edu/J ACpublic/JCMT/ 



ISchoier et al.l d2006bl) and iRamstedt et all d2006l 120081 
The abundance distribution o f CO is based on the photochem- 
ical models of iMamon et al.l d!988l) adopting an initial (pho- 
tospheric) abundance of 1 x 10~ 3 for carbon stars, 6 x 10~ 4 for 
S-type AGB stars and 2 x 10~ 4 for M-type AGB stars. The ki- 
netic temperature structure is obtained in a self-consistent man- 
ner from solving the energy balance equation, where the CO line 
cooling is directly obtained from the excitation analysis. Mass- 
loss rates derived from detailed radiative transfer modelling of 
CO multi-transition line e mission are accurate to within a factor 
2-\ dRamstedt et alj|2008l) . 

The excitation analysis includes the first 41 rotational lev- 
els within the ground state and the first vibrationally excited 
(v= 1) state of CO at 4.6 yum. The collisional rates are taken 
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treated as described in Lindqvist et al. (2000). Frequencies for 
the hyperfine transitions are obtained f rom the CDMS database 
dMtiller et al.ll200lHMuller et al.ll2005l) . 

When the molecule is bending and rotating simultaneously, 
the double degeneracy of the bending mode will be lifted, re- 
sulting in the splitting of the bending mode (010) into two lev- 
els (01 Ic and Ql ld Q). This is referred to as /-type doubling 
and is included in the excitation analysis. The transition prob- 
abilit ies and frequencies are taken from the GEISA database 
dJacquinet-Husson et al . 2001. 

The collision al rate coefficients are taken from 
iGreen & Thaddeusl d 19741) and G reen (unpublished dat£0). 
The rate coefficients presented in IGreen & Thaddeusl ([l974) 
were calculated for temperatures in the range from 5 to 100K, 
including energy levels up to 7 = 7 for collisions with He. The 
unpublished rates by Green were obtained for the lowest 30 
rotational transitions and for temperatures from 100 to 1200 K, 
again for collisions with He. These two sets of rate coefficients 
were subsequently scaled by a factor 1.37 to represent collisions 
with H2. Extrapolation to include all energy levels up to J = 29 
for temperatures between 5 and 120 0K was made and can be 
obtained from the LAMDA database (Schoi er et al.l2005l). Quite 
recent ly, new collision rates were published bv lDumouchel et~al1 
J20T0I) . We tested using these rates for three objects with low, 
intermediate, and high mass-loss rate. The results differed by 
less than a few percent compared t o the results when u s ing th e 
scaled and extrapolated data from IGreen & Thaddeusl dl97l . 
for all three cases. 

The results from the HCN excitation analysis are given in 
Tables[3]and0 



4. Results and discussion 

4.1. Results from the radiative transfer modelling 

In order to determine accurate abundances it is important to 
know the spatial extent of the molecular envelope. Since dif- 
ferent molecular transitions have different excitation conditions, 
multi-transition observations have the potential to constrain both 
the initial abundance, fo, and e-folding radius, r e (Eq. |2]i and 
hence describe the distribution of HCN molecules in the CSE. 
The best-fit initial abundances provide constraints on chemical 
models in and near the photosphere as well as dust condensa- 
tion models. The envelope sizes obtained, on the other hand, 
co nstrain photochemical m o dels o f the CSE. Thi s is illus trated 
Gonzalez Deleado et all d2003h . ISchoier et all (l2006bl) and 



Ramstedt eTalJ d2009l) for SiO. and in ISchoier et al J <l2007al) for 
SiS. 

Due to lack of observational constraints a full analysis, deter- 
mining both the initial abundance and the envelope size, is possi- 
ble only for a sub-sample (20) of sources. For these sources, both 
parameters are considered as free parameters and varied within 
reasonable ranges until the best fit to the observations are found. 
The results are then used to derive a relation describing how the 
HCN envelope size scales with the wind density (M/v e ). This re- 
lation is used to estimate the HCN envelope size for all sources 
and the initial abundance is the only free parameter in the ra- 
diative transfer analysis. Given the uncertainties involved in the 
modelling of carbon stars as explained further in Sect. I4.1.2l (par- 
ticularly for the J -I — >0 transition) no carbon stars are used 
when deriving the HCN envelope size relation. 



from lFlowerl (l200ll) and IWernli etail (120061) . They were calcu- 
lated for temperatures in the range from 5 to 400 K including 
rotational energy levels up to J — 29 and J = 20 for collisions 
with para-H2 and ortho-H2, respectively. The collisional rate co- 
efficients have been extrapolated to include energy levels up to 
J = 40 and collisional temperatures up to 2000 K, as described 
in lSchoier et al.1 (120051) . and are made publicly available through 
the Leiden Atomic and Molecular Database (LAMDAfl An 
ortho-to-para ratio of 3 was adopted when weighting together 
collisional rate coefficients for CO in collisions with ortho-H2 
and para-H2. 

The parameters obtained from the CO excitation analysis are 
reported in Tables |3]and|4] 

3.4. HCN line modelling 

Once the basic envelope parameters are known from the CO line 
modelling and dust radiative transfer, the abundance of HCN 
in the circumstellar envelopes can be determined. The HCN 
abundance distribution (i.e. the ratio of the number density of 
HCN molecules to that of H2 molecules as a function of radius), 
/ = n(HCN)/n(H2), is assumed to be described by a Gaussian 







2\ 


- 










t 



where /o is the initial abundance, and r e the e-folding radius for 
HCN. ' 

HCN is a linear molecule and has three different vibrational 
modes: the CN stretching mode vi, the CH stretching mode V3, 
and the doubly degenerate bending mode V2 (doubly degenerate 
since the bending of a non-rotating linear molecule is direction 
independent and will result in radiation of the same frequency). 
The excitation analysis includes radiative excitation through the 
CH stretching mode (v3 = l; also denoted 001) at 3yt/m, and 
the bending mode (V2 = 1; also denoted 010) at 14ytzm. The CN 
stretching mode (vi = 1; also denoted 100) at 4.8 //m has transi- 
tions that are about 20-1000 times weaker and it is therefore 
not included in the analysis. In each of the vibrational levels we 
include rotational levels up to J — 29. 

The first vibrational overtone (020) lies at energies above 
the ground state corresponding to 7 /urn and has transitions to 
the (010) and (000) states with Einstein A coefficients of the 
same order as those from the (010) to the (000) state (e.g., 
IZiurvs & Turner1ll986l) . A test was performed to estimate its in- 
fluence on the excitation analysis covering optical depths that 
are characteristic for the sample sources. For all models, the to- 
tal integrated intensities of the lines included in this paper are 
changed by less than a few percent when the (020) state is in- 
cluded. At high optical depths the radiative excitation is of mi- 
nor significance and therefore likely also the excitation to higher 
vibrational states. An effect on higher rotational lines at interme- 
diate and low optical depths, cannot be excluded, but to evaluate 
this is beyond the scope of this paper. For the excitation anal- 
ysis and intensities of the lines used to constrain the modelling 
performed in this paper, we conclude that excitation through the 
(020) state will have no significant effect. 

The nuclear spin of the 14 N atom results in the splitting of 
the rotational levels into three hyperfine components. Hyperfine 
splitting is included only for rotational levels in the ground state 
up to J — 4 (where the splitting is of the same order as the lo- 
cal line width) and the resulting line overlaps are accurately 

3 http : //www . strw . leidenuniv . nl/ ~moldata 



4 http : //data . giss . nasa . gov/mcrates/data/hcn_he_rates . txt 
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Table 2. HCN radial abundance distribution results. 
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= 3.4 


X 


io- 


•6 


2.2 = 


= 1.7 


X 


10 16 


3 


4.5 


M-type 
























RXBoo 


4.7 = 


= 2.3 


X 


io- 


-7 


3.5 = 


= 1.4 


X 


10 15 


6 


2.2 


TXCam 


2.8 = 


t 1.1 


X 


io- 


-7 


4.5 = 


= 2.7 


X 


10 16 


7 


0.8 


RCas 


3.5 = 


t 1.4 


X 


io- 


-7 


5.9 = 


= 2.1 


X 


10 15 


6 


3.5 


R Dor 


9.0 = 


= 4.0 


X 


io- 


-8 


7.1 = 


= 4.0 


X 


10 15 


3 


2.0 


WHya 


8.5 = 


= 4.5 


X 


io- 


-7 


3.0 = 


= 1.1 


X 


10 15 


4 


6.1 


RLeo 


2.5 ± 1.1 


X 


io- 


-7 


6.0 = 


= 2.5 


X 


10 15 


4 


0.9 


WXPsc 


3.0 = 


t2.0 


X 


io- 


-7 


3.0 = 


= 2.0 


X 


10 16 


5 


1.0 


IKTau 


4.3 = 


t 2.8 


X 


io- 


-7 


1.4 = 


= 0.8 


X 


10 16 


4 


0.5 


IRC- 10529 


4.5 = 


= 2.5 


X 


io- 


-8 


3.2 = 


t 1.8 


X 


10 16 


4 


2.7 


IRC-30398 


1.1 = 


= 0.3 


X 


io- 


7 


2.6 = 


t 1.3 


X 


10 16 


4 


0.9 



" The abundance distribution is assumed to be Gaussian (see Eq.0 



4.1 .1 . HCN envelope size 

In the present sample 12 M- and 8 S-type AGB stars have been 
observed in three or more different lines. For these sources, the 
size of the HCN radio-line-emitting region (r e ) and the inner 
HCN abundance (fo) are varied simultaneously. In 16 cases it 
is possible to constrain both parameters reasonably well, as re- 
ported in Table [2] The sensitivity of the line emission to vari- 
ations of the two adjustable parameters fo and r e is illustrated 
in Fig.[TJfor the high-mass-loss-rate M-type AGB star TX Cam 
and the intermediate-mass-loss-rate S-type AGB star W Aql. In 
Fig- El the model spectra are directly compared with the obser- 
vations for the best-fit models. The results for the intermediate- 
mass-loss-rate carbon star R For are also shown in Figs [TJand [7] 
as an example, however they are not included when deriving the 
HCN envelope size relation. In most cases the reduced^ 2 values 
ttred = XioiKN - 2)) are < 3 indicating reasonably good fits to 
the observed HCN line emission. 

In Fig. [2] the fractional abundances for this sub-sample are 
plotted as a function of a density measure (M/v e ). No apparent 
trends with density are present. However, the S-type AGB stars 
appear on average to have larger HCN fractional abundances 
than M-type AGB stars (this is further discussed in Sect. 14.1.31 ). 
The carbon star R For is again shown for comparison. 

We find that the derived envelope sizes in Table[2]scale with 
density as 



M 

logr e = 19.9(± 0.6) + 0.55(± 0.09) log| — |, 



(3) 



where v e is the expansion velocity of the wind in km s , and M 
the mass-loss rate in M yr _1 (the values in parenthesis are 1 cr 
errors). As shown in Fig. [2] the scaling law gives a good estimate 
of the envelope sizes, within the uncertainties of the results from 
the multi-transition radiativ e transfer. The reduced y 2 of th e fit 
is 1.5 indicating a good fit. iGonzalez Delgado et al.l d2003b de- 
rived a slope of 0.48 in their analysis of SiO line emission which 



^10 




10 10 

M/v e [Mq yr" 



10 

1 km" 1 



10 



-5 



Fig. 2. Results from excitation analysis of HCN multi-transition 
single-dish observations assuming a Gaussian abundance distri- 
bution (see Eq. Top panel - Derived initial fractional abun- 
dance (fo) as a function of envelope density (M/v e ) for M- 
(squares) and S-type AGB stars (triangles). Bottom panel - 
Derived envelope sizes (r e ) as a function of envelope density 
(M/v e ) for M- (squares) and S-type AGB stars (triangles). The 
solid line is a least-squares fit to the data points as described by 
Eq.[3] The dashed and dotted lines show the results from a simple 
photochemical model (Sect. l4~7TTTT i. The dashed line is calculated 
for h-0.2 and the dotted line for /i=0.5 (see text for details). The 
carbon star R For (dot) is a also shown in the plots. 



is consistent with the value of 0.55 obtained here, within the un- 
certainties. The HCN sizes derived here are systematically larger 
than the corresponding SiO envelope sizes, although they are 
still comparable within the mutual uncertainties of the scaling 
laws. 

High angular resolution HCN observations, capable of re- 
solving the emitting regions, require the use of interferometers 
and exist only for a small number of sources. Such observations 
generally put better constraints on the HCN abundance distribu- 
tion than multi-transition single-dish observations and in order 
to test our relation (Eq. |3) w e have compar ed it to previous in- 
terferometer results. Lindqvist et al. (2000) estimated envelope 
sizes for a sample of five carbon stars with varying mass-loss 
rates. In three cases it was possible to obtain good constraints. 
For the carbon stars CW Leo and LP And, Eq. [3] gives estimate s 
that agree very well with the findings of Lindqvi st et alJ (feOOO). 
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Fig. l.^ 2 -map showing the quality of the fit to available HCN multi-transition single-dish data for the AGB stars R For (left panel), 
W Aql (middle panel), and TX Cam (right panel) when varying the adjustable parameters, /o and r e , in the model. Contours are 
drawn at the 1, 2, and 3 cr levels. Indicated in the upper right corner is the reduced x 2 of the best-fit model and the number of 
observational constraints used, N. 



In the case of RW LMi, Eq.[3]gives an envelope size that is about 
a factor of two too small. 

In Fig. [3] we show the interferometric observations of HCN 
J = 1 — > emission toward s the two M-typ e AGB stars TX Cam 
and IK Tau, performed by iMarvell (120051) at the Owens Valley 
Radio Observatory (OVRO). Plotted with the data are model re- 
sults for different HCN abundance distributions within reason- 
able ranges of the results found in Table [2] The analysis and 
comparison is performed in the Mv-plane in order to maximize 
the sensitivity and resolution of the observations. It is found that 
envelope sizes of 1.6 x 10 16 cm for TX Cam and 9.0 x 10 15 cm 
for IK Tau provide good fits to the observations. These values are 
consistent with those obtained from the multi-transition analysis 
of single-dish data reported in Table [2] Furthermore, a compar- 
ison between the line strength of the HCN / = 1 — » emission 
measured with the OSO 20 m telescope, and the sum of the emis- 
sion measured by the OVRO interferometer, shows good agree- 
ment, and no indication of extended emission being resolved out 
by the interfe rometer. 

Recently, iMuller et al.1 (120081) reported interferometric ob- 
servations of HCN ./ = 8 — > 7 emission towards the low-mass- 
loss-rate M-type AGB star W Hya. They find, from a de- 
tailed radiative transfer analysis, an HCN envelope size of about 
3 x 10 15 cm, in very good agreement with our derived value of 
3.0+ l.lx 10 15 cm (see Table |2j. These examples give us con- 
fidence in the adopted approach and we feel that the envelope 
sizes, as estimated by Eq. [3] are adequate for use in a statistical 
study of a large sample of sources such as this. 

The HCN envelope sizes can also be estimated using a sim- 
ple photochemical model. Assuming that the siz e is regulated by 
photodissociation (Hug gins & Glassgoldlll982h . the abundance 
distribution is given by: 

df = -G exp(-?) 
dr v exp 

where Go is the unshielded photodissociation rate of HCN 
and rd is the dust-shielding distance: 

. „ 3 gabs M d 

r d = 1.4- . (5) 

4adp d 4nv d 

<2abs is the dust absorption efficiency, p d and a d are the dust 
grain density and radius, and M d and v d are the dust mass- 
loss rate and veloc i ty. By using the h-parameter defined in 
ISchoier & Olofssonl (120011) . and assuming ¥ = 0.01, p d = 



2 g cm 3 , Ad = 0.05 /jm, the dust-shielding distance can be re- 
written as: 



r d = 5.27 10 22 cm, (6) 

Vd 

The wavelength dependence of <2abs in the region of interest, 
-1000-3000 A, is weak (e.g.. ISuhl l2000l) and a generic value 
of 1.0 is adopted. As a comparison to our observational results, 
the e-folding radius calculated from Eq. |4] is shown in Fig. [2] 
for h = 0.2 (dashed) and h = 0.5 (dotted) (typical values for low- 
, L < 5000 L , and high-luminosity, L > 500 L , stars) for 
different M/v e . Go is taken to be l.lx 10 _9 s _1 (Ivan Dishoeckl 
1988), and the dust velocity is calculated from the gas veloc- 
ity, measured from the width of the CO lines, and the drift ve- 
locity, calculated from the luminosity and the mass-loss rate 
( ISchoier & Olo fsson 2001). The results from the photochemical 
model is further discussed in Sect. 14.21 

4.1.2. HCN maser action 

The HCN J = 1 — > line has been found t o exhibit maser fea- 
tures in a number of low-mass carbon stars (Olof sson et al.l!993l 
Il998h . ILindqvist et al.1 (l2000h found that their radiative transfer 
model could not properly account for the observed intensities in 
this type of sources. The observed intensities were significantly 
higher than tho se predicted by their model. T his was further cor- 
roborated by iDinh-V-Trung & Rieul (l2000|) in their modelling 
of the low-mass-loss-rate object Y CVn. ILindqvist etaO {2000) 
speculated that maser action, emphasized by deviations from a 
smooth wind, could be the re ason for the discrepanc y from the 
model predictions. Recently, Shinnaga et al. (2009) have also 
found evidence for HCN maser emission in the inner wind of 
the high-mass-loss rate object CW Leo (IRC +10216), suggest- 
ing that maser emission is a common phenomenon in carbon 
stars. 

In order to investigate where this phenomenon arises, we 
have calculated a number of models covering the density range 
of our sample sources. In these test models the mass-loss rate 
is varied in the range from 1 x 10" 7 M yr _1 to 3 x 10~ 5 M yr" 1 
for a wind that expands at a velocity of lOkms 1 . The central 
star is assumed to have a temperature of 2400 K and a luminos- 
ity of 5000 L in all models. The CO abundance is assumed to 
be 6 x 10~ 4 . To cover the majority of the stars in our sample, two 
different HCN abundances; 1 x 10~ 5 and 1 x 10 -7 , were consid- 
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equivalent spatial scale ["] 




equivalent spatial scale ["] 
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Fig. 3. Visibilities, averaged over 2km s _1 around the sys- 
temic velocity, obtained for IK Tau and TX Cam using 
the Owens Valley Millimeter Array (iMarvell 120051) . The 
observations are overlaid by model results using various 
abundance distributions. Top panel - For TX Cam the follow- 
ing models are shown: /o=l.lxlO~ 6 and r e = 1.2 x 10 16 cm 
(dotted line), / = 7.0xl(T 7 and r e = 1.6 x 10 16 cm (solid 
line), and / = 5.0xl(T 7 and r e = 2.0 x 10 16 cm (dashed 
line). Bottom panel - For IK Tau the following models 
are shown: /o = 1.0xlO~ 6 and r e =6.0x 10 15 cm (dotted 
line), /o = 5.0xl(T 7 and r e = 9.0 x 10 15 cm (solid line), and 
fo = 3.0 x 10~ 7 and r e = 1.2 X 10 16 cm (dashed line). 



ered. Under certain assumptions, Eq. [4]c an be solved to give an 
estimate of the photodissociation radius (Olofss on et alJ l998): 



r ph = 10 s l^-l cm i 2.9 x 10' 



M 



lio) ' 



cm, 



(7) 



where v e is the expansion velocity of the wind in kms -1 , and 
M the mass-loss rate in M yr _I . This estimate of the enve- 
lope size was then further scale d b y a factor of two, as sug- 
gested by Lindq vis t et all d2000t) and lBieging et al.1 12000) (see 
also Sect. 14. 1 .Ti l. The source was placed at a distance of 500 pc 
and the beams for the HCN lines were assumed to conespond to 
those from a 15 m telescope, such as the JCMT and the SEST. 
All observational results shown in Fig. [4] have been scaled to a 



Z 
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Fig. 4. Observed (markers) and modelled (lines) intensity ratios 
between the HCN J = 1 — > and J = 3 — > 2 lines as a function 
of a density measure (M/v e ) of the CSE. Observations of car- 
bon stars are indicated by blue circles, M-type AGB stars by 
red squares, and S-type AGB stars by green triangles. Model 
predictions with HCN abundances of 1 x 10~ 5 (solid line) and 
1 x 10~ 7 (dashed line) are shown for comparison for a circum- 
stellar model as described in the text. 



15 m telescope assuming unresolved emission (i.e., that the in- 
tensity scales as the inverse of the aperture squared). The uncer- 
tainty in the observed intensity ratios is estimated to be « 50%. 

The model line intensity ratios between the J -3—>2 and 
J = 1 — » lines (dashed and dotted lines, one for each assumed 
HCN abundance) are compared to observed values (markers) 
in Fig. |4] There appears to be a distinction between the car- 
bon stars and M-type AGB stars in that carbon stars tend to 
have lower HCN J -3—>2 to 7=1— >0 line intensity ratios in 
the low-density regime. The S-type AGB stars tend to fall in 
between the two other chemical types. We identify a region 
where M/v e < 3 x 10~ 8 M yr" 1 km" 1 s and the HCN J = 3 -» 2 
to J — 1 — » line intensity ratio is ~ 1 for several carbon stars. 
Within the basic assumptions of the circumstellar model it is 
not possible to achieve such a low line intensity ratio. We be- 
lieve these low values are due to maser action (possibly strength- 
ened by a clumpy gas structure) in the J = 1 — > transition, con- 
cluding that under the assumptions of our circumstellar model it 
will be very difficult to reproduce the observations of the lower- 
density carbon stars where the conditions for strong maser action 
are met (high abundance/low density). In the higher density case, 
the maser action is likely more effectively quenched and there- 
fore the higher-density carbon stars are not subject to the same 
problems. When comparing model predictions with observations 
in this case one should keep in mind that a relatively large scat- 
ter is expected around the lines representing model results, due 
to significant deviations in individual sources from these model 
assumptions. Therefore, Fig. |4]is not giving strict limits for the 
radiative transfer model, but illustrates the difficulties and em- 
phasizes that specific care has to be taken when modelling and 
interpreting the results on carbon stars in this region of parame- 
ter space. In the following analysis, the modelling of these stars 
is performed excluding the J =1 — > transition when the (bad) 
fit to it, dominates the totals-estimate. 
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4.1.3. HCN abundances 

Using the scaling relation in Eq. [3]it is now possible to estimate 
the HCN envelope sizes for all of the sample sources. Once the 
size, r e , is known only one free parameter remains in the excita- 
tion analysis, the initial abundance fy. The derived abundances 
and adopted envelope sizes of HCN are reported in Tables [3] 
and 3] and Figs [5] and [6] The fits to the data are generally good 
with reduced^ 2 values of « 1-3. It is found that the abundance 
of HCN varies substantially with the photospheric C/O-ratio of 
the central star. We find that the median HCN abundance is 
2.9 x 10~ 5 for the carbon stars (25 stars; dash-dotted line in Fig.|5] 
and filled circles in Fig.[6]l, more than two orders of magnitudes 
higher than for the M-type AGB stars (median is 1.2 x 1CT 7 ; 25 
stars; dashed line in Fig. [5] and filled squares in Fig. [6). The S- 
type AGB stars, on the other hand, seem to be distributed in be- 
tween the two other types (19 stars; solid line in Fig.|5]and filled 
triangles in Fig. [6j with a median HCN abundance of 7.0 x 10~ 7 . 
The spread in abundances (as measured by the ratio of the 75 th 
and 25 th percentile) is also larger for the S-type AGB stars (5.6) 
compared to those of the carbon stars (4) and the M-type AGB 
stars (2), possibly indicating the sensitivity of the HCN abun- 
dance in the border area between an oxygen-dominated chem- 
istry (lower HCN abundances) and a carbon-dominated chem- 
istry (higher HCN abundances) within the S-type sample. 
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Fig. 5. Histograms showing the derived HCN fractional abun- 
dances (/o) for the carbon stars (dash-dotted, blue line; 25 stars), 
M-type AGB stars (dashed, red line; 25 stars), and S-type AGB 
stars (solid, green line; 19 stars). 

There is no apparent trend that the HCN abundance de- 
creases as the density of the wind (M/v e ) increases, as would 
be expected if condensation of HCN molecules onto dust grains 
were an important process in the dense inner parts of the wind. 
The Pearson correlation coefficients are r — - 0.03 for the carbon 
stars, r = 0.49 for the M-type AGB stars and r = - 0.36 for the S- 
type AGB s tars. This differs from the result of our previous anal- 
ysis of SiO dGo nzalez D elgado et alJl2003HSchoier et alj2006bt 
Ramstedt et alj|2009l) w here clear cor r elation s with the density 
of the wind are present. ISchoier e t al. (2007a|) also found weak 
correlations of SiS fractional abundances with M/v e for a sample 
of carbon stars and M-type AGB stars. 

The uncertainty in the derived abundances is about a factor 
of 3 within the adopted model (see Table 0. Additional uncer- 
tainties are introduced by the simplifications of the circumstellar 
model. A systematic error due to the adopted CO/H2 abundance 



might also be p resent, but we e stimate this to be no larger than a 
factor of -1.5 (Ramsted t et alJl2008l) . The estimated uncertain- 
ties cannot account for the differences that we find between the 
chemical types. Also, we believe that the spread among the S- 
type stars is real and indicative of the dependence of the HCN 
abundance on the C/O-ratio of the star. 

The HCN abundances that we obtain ag ree quite well , within 
a factor of four, with those presented by IWoods ef al. (2003) 
using a simple excitation model, for the five high-mass-loss- 
rate carbon stars in common. The abundances derived in the 
present analy sis are usually higher than the ones reported by 
I Woods et a l. (2003). This discrepancy may be a result of their 
adopted excita tion temperature being too high and of optical 
depth effects. Buiarraba fet al.1 (119941) found HCN abundances 
that on average are ~ 20 times higher for the carbon stars than 
for the M-type AGB stars, support ed by the finding s in th is pa- 
per. The average abundances that iBuiarrabal et al.l (11994b de- 
rived, through a simple analysis, are 7.4 x 10 -6 for carbon stars, 
3.7 x 10~ 7 for M-type AGB stars and 4 x 10~ 6 (only three 
sources) for S-type AGB stars. IZiurvs et al.l(l2009l) derived HCN 
abundances for five oxygen-rich stars using new NRAO data, 
together with literature data and a model utilizing the Sobolev 
approximation. For W Hya the results agree very well, but for 
TX Cam and IK Tau, the differences are too big to be easily 
explained. Since both sources are high-mass-loss-rate sources it 
is po ssible that the discrep ancy is due to optical depth effects. 
For iLindqvist et al.l (l2000h modelled, spatially resolved, inter- 
ferometric observations of HCN ./ = 1 — > line emission from 
a sample of five carbon stars and derived abundances that all 
agree within a factor of two with our results for the correspond- 
ing sources. 

4.2. Comparison with chemical models; HCN abundances 

The abundance of HCN derived using a chemical model 
assuming thermal equilibrium (TE) is highly dependent on 
the assumed C/O-ratio of the gas. For a typic al M-type 
star (IK Tau with C/O = 0.75, iDuari et al.l (Il999h ) the ex- 
pected TE abundance is more than six orders of magnitude 
smaller than th at derived for a typical carbon star [IRC +102 16 
with C/O = 1.5. IWillacv & Cherchnefij (Il998l) l. For S-type stars 
(0.98 <C/0< 1.01) the HCN abundance is expected to have 
a spread of about three orders of magnitude assuming TE 
(ICherchneffll2006l) . 

The inclusion of shocks in the chemical models drastically 
alters the expected abundance of HCN, particularly for the 
oxygen-rich case (Cherchneff 2006|). For a C/O-ratio of 0.75, 
the abundance is changed from 2xl0" n to lxlO 5 at 5 R* when 
shocks are included. For C/0=1.5, the HCN abundance is ex- 
pected to be lower (3xl0~ 6 ) in the shock model than both the 
TE value (5xl0~ 5 ) and that derived for the M-type stars. For the 
S-type stars the spread in the abundances when varying the C/O- 
ratio around 1 is less than a factor of two, and for C/0=1 an 
abundance of about 6.5xl0~ 6 is expected at 5 R*. 

The derived circumstellar abundances in this work are, on 
average, slightly less than four orders of magnitude larger for 
the M-type stars than the expectations from TE models, while 
the results for the carbon stars agree quite well with TE expec- 
tations. The spread in derived abundances around the average 
is small for both types of stars. Compared to the results from 
shock chemistry models, our derived abundances are on aver- 
age about an order of magnitude larger for the carbon stars. A 
discrepancy of the same order has previ ously been noted for 
IRC+10216 dWillacv & Cherchneffl[l998l) . For the M-type stars 
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Table 3. Model results. 



SED modelling CO modelling HCN modelling 
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Fig. 6. Results from multi-transition excitation analysis of circumstellar line emission towards a large sample of AGB stars with 
varying physical and chemical properties, left - HCN fractional abundance, /o(HCN), as a function of a density measure (M/v e ) 
for carbon stars (filled circles), M-type AGB stars (filled squares) and S-type AGB stars (fille d triangles), middle - SiO fractional 
abundance, /o(SiO), a s a function of a density measu re (M/v e ) for carbon stars (filled circles; Schoi er et al.lE0 06b). M-type AGB 
stars (filled squares; iGonzalez Delgado et al.l 12003) and S-type AGB stars (filled triangles; Ramstedt et al. 2009). right - SiS 
fractional abu ndance, fo(SiS), as a function of a density measure for carbon stars (filled circles) and M-type AGB stars (filled 
squares) from Sch oier et al.l (l2007al) . 



our results are, on average, approximately one order of magni- 
tude le ss than the theore tical value from the models including 
shocks dCherchneffl2 006 ) . For the S-type stars, the derived abun- 
dances are, on average, about an order of magnitude larger than 
the TE value. Notably, for these stars, the spread in circumstellar 
abundances is relatively large. As discussed below (Sect. 14.31 ), 



abundances derived from chemical models including shocks are 
sensitive to the assumed sho ck parameters, i.e., the shock ve- 
locity (e.g. iDuari et al.| [l999). Since the observed sources are 
expected to cover a large range in shock parameters, while the 
models are often run for one set of parameters, this can some- 
times limit the comparison. However, the clear dependence of 
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Table 4. Model results. 
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the circumstellar HCN abundances on the C/O-ratio is in contrast 
with the expectations from chemical models including shocks, 
where a homogeneous chemistry is derived with respect to the 
C/O-ratio. Also, the observational HCN abundances show no, 
or a very weak, dependence on the wind density, again indicat- 
ing that the chemistry is not influenced by shocks. Our results 
are also in contrast with the conclusions of iDecin et al.l ([2008) 
where chemical homogeneity (regardless of C/O-ratio and evo- 
lutionary stage) is deduced from observational data (two carbon 
stars, one S-type star, and one M-type star). 



The reason for the discrepancy between the chemical mod- 
els and the observational results is difficult to identify with any 
certainty. Given the high abundances found for the M-type stars, 
non-equilibrium chemical processes (such as shocks) certainly 
seem to play a role in the formation of HCN. A common for- 
mation route for HCN is through the react ion between CN and 
H 2 forming HCN and H (iDuari et al.lll999h . The rate for this re- 
action is highly temperature dependent and therefore the abun- 
dance of HCN found in the chemical models will be very sen- 
sitive to the assumed shock properties. However, this would 
not explain the differences between our results for the differ- 
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ent chemistries and the resul ts from chemica l models when the 
same shock model is used (ICherchn eff 2006). Given that CN is 
important for the formation of both HCN and CS, a similar in- 
vestigation as the one presented here, of the abundances of CN 
and CS would be very helpful for the further interpretation of 
our results. A complicating factor could be, however, that in the 
CSEs, CN is also produced through photodissociation of HCN. 
Previous observations indicate that the abundance of CS is also 
highly dependent on the C/O-ratio (lOlofsson et al .1119981) . 



4.3. Observational circumstellar abundances 

In Fig. [6] we present the HCN results together with the 
results o f similar abundance analysis on the sam e samples 
for SiO dSchoier et all l2006bt iRamstedt etail |2009), and SiS 
dSchoier et al.ll2007ah . The results are shown as functions of the 
circumstellar density. The SiS analysis has only been performed 
on a sub-sample of carbon and M-type stars. 

It is clear that the formation of both HCN and SiS are depen- 
dent on the C/O-ratio of the star, although less than what would 
be expected from TE models. The HCN results are also consis- 
tent with the the S-stars constituting a distinct sample chemi- 
cally. From the theoretical models, SiS is expected to have a de- 
pendence on the C/O-ratio, also when chemical non-equilibrium 
processes are taken into account. SiO, on the other hand, shows 
no apparent dependence on the C/O-ratio, neither observation- 
ally, nor theoretically. According to current chemical models, 
HCN and SiO are expected to show the same behaviour, and the 
reason why HCN is significantly more sensitive to the C/O-ratio 
is not clear. 

In a shock-induced chemistry, the resulting abundances of 
all molec ules are expected to be sensitive to the assumed shock 
velocity (Cherchneff 2012). Although the exact dependence is 
not entirely straight-forward, a correlation with the pulsation pe- 
riod of the star might therefore be expected. However, no such 
dependences are found for HCN and SiS, and the one found for 
SiO probably reflects the correlation between mass-loss rate and 
pulsation period. This attempt to correlate the estimated abun- 
dances with the shock properties is rather crude, because any 
radial dependence will be smeared out. There are also obser- 
vational indications that the abundance distribution of all three 
molecules is more complicated than assumed in our radiative 
transfer m odel for many sources (see e . g., lin e-profile shapes in 
Fig. 171 and lSchoier et al.ll2006ai l2007al (20111) . 

The theoretically predicted abundances do not take conden- 
sation onto grains into account, while the observational results 
apply to the post-condensation region. Hence a comparison be- 
tween them may yield some insight into dust grain formation and 
evolution. For SiO there is a clear dependence on the circumstel- 
lar density and this is interpreted as due to SiO adsorption onto 
dust in the high-density wind. On the contrary, the HCN and the 
SiS abundances show no, or only very weak, dependence on the 
circumstellar density, in the HCN case over more than two or- 
ders of magnitude in density. From this one may conclude that 
these species are not affected by grain formation and/or evolu- 
tion. However, we caution that a comparison between the theo- 
retical and observational values is hampered by the significant 
uncertainties in both models, and by the fact that the chemical 
theoretical models are usually calculated for only one set of stel- 
lar parameter, while the observed samples cover large ranges of 
stellar mass, luminosity, periods, shock properties, and tempera- 
tures. 



5. Conclusions 

New multi-transition (sub-)millimetre HCN line observations of 
a statistically significant sample of AGB stars with varying pho- 
tospheric C/O-ratios and circumstellar properties are presented. 
From a detailed excitation analysis, based on a reliable physical 
model of the sources, we reach the following conclusions: 

- The median derived abundances of HCN are 2.5 x 10 5 , 
7.0 x 10~ 7 , and 1.2 x 10~ 7 for the carbon stars (25 stars), S- 
type AGB stars (19 stars) and M-type AGB stars (25 stars), 
respectively. 

- We derive a scaling law for the size of the HCN molec- 
ular envelope as a function of the wind density (M/v e ). 
The envelope sizes estimated using this law are found to 
agree with previous interferometric observations and with 
a simple photochemical model. The HCN envelopes are 
systematically larger than the corresponding SiO envelopes 
(Gonzalez Del gado et al.l2003h . however still within the mu- 
tual uncertainties of the scaling laws. 

- There is a clear dependence of the derived circumstellar 
HCN abundances on the C/O-ratio of the star. Carbon stars 
have, on average, about two orders of magnitude higher 
abundances than M-type AGB stars. The S-type AGB stars, 
which appear to have a larger spread in their derived HCN 
abundances, typically fall in between the two other types, 
however, slightly biased towards the M-type AGB stars. 

- The sensitivity of the circumstellar HCN abundance to the 
photospheric C/O-ratio is in star k contrast to pre dictions 
from non-LTE chemical models (ICherchneffl 120061) . where 
very little difference is expected in the HCN fractional abun- 
dances between the various types of AGB stars. 
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Fig. 7. Multi-transition spectra (histograms) of HCN line emission (HHT data from Bieging et al. 2000), overlaid by spectra from 
the best-fit model (solid lines; assuming a Guassian abundance distribution), towards top panels - the carbon star R For using a 
fractional HCN abundance /o = 4.0 x 10~ 5 and an envelope size r e = 9.3 x 10 15 cm. middle panels - S-type AGB star W Aql using 
a fractional HCN abundance /o = 1.2 X 10" 6 and an envelope size r e = 1.3 X 10 16 cm. bottom panels - M-type AGB star TX Cam 
using a fractional HCN abundance /o = 3.5 x 10~ 7 and an envelope size r e = 2.0 x 10 15 cm are also shown (solid lines). The hyperfine 
splitting of the J — 1 — > transition significantly broadens the line. This effect is explicitly taken into account in the modelling. The 
calibration uncertainty in the observed spectra is » ± 20%. 
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Appendix A: New observations of HCN line 
emission 

The spectra obtained from our new multi-transition HCN survey 
are shown in Figs IA.114A.4l See Sect. [2] for details on the ob- 
servations. The corresponding (velocity) integrated intensities, 
J r m b dv, are reported in Table lA.fl & IA. 21 including also previ- 
ously detected lines that are used in the HCN excitation analysis. 
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Table A.l. Integrated (7 b s = J r m b dv) line intensities in K km s _1 for the HCN (v = 0, / — » /— 1) line emission used in the excitation 
analysis. 
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Table A.2. Integrated (7 b s = J r m b dv) line intensities in K km s _1 for the HCN (v = 0, / — » /— 1) line emission used in the excitation 
analysis. 
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A.l. New multi-transition HCN observations of carbon stars. 
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. A.2. New multi-transition HCN observations of carbon stars. 
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Fig. A.3. New multi-transition HCN observations of S-type stars. 
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